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Abstract 

Magnetic resonant coupling (MRC) is a practically appealing method for realizing the near-held wireless power 
transfer (WPT). The MRC-WPT system with a single pair of transmitter and receiver has been extensively studied 
in the literature, while there is limited work on the general setup with multiple transmitters and/or receivers. In this 
paper, we consider a point-to-multipoint MRC-WPT system with one transmitter sending power wirelessly to a set 
of distributed receivers simultaneously. We derive the power delivered to the load of each receiver in closed-form 
expression, and reveal a “near-far” fairness issue in multiuser power transmission due to users’ distance-dependent 
mutual inductances with the transmitter. We also show that by designing the receivers’ load resistances, the near-far 
issue can be optimally solved. Specifically, we propose a centralized algorithm to jointly optimize the load resistances 
to minimize the power drawn from the energy source at the transmitter under given power requirements for the 
loads. We also devise a distributed algorithm for the receivers to adjust their load resistances iteratively, for ease of 
practical implementation. 


Index Terms 

Wireless power transfer, magnetic resonant coupling, multiuser charging control, optimization, iterative algo¬ 
rithm. 


I. Introduction 

Inductive coupling iCQ-El is a conventional method to realize the near-field wireless power transfer (WPT) 
for short-range applications up to a couple of centimeters. Recently, magnetic resonant coupling (MRC) If5l-f7t 
has drawn significant interests for implementing the near-field WPT due to its high power transfer efficiency for 
applications requiring longer distances, say, tens of centimeters to several meters. The transmitter and the receiver 
in an MRC-WPT system are designed to have the same natural frequency as the system’s operating frequency, 
thereby greatly reducing the total reactive power consumption in the system and achieving high power transfer 
efficiency over long distances. 

The MRC-WPT system with a single pair of transmitter and receiver has been extensively studied in the literature 
for e.g. maximizing the end-to-end power transfer efficiency or the power delivered to the receiver with a given 
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input power constraint ||8l- llTTl . However, there is limited work on analyzing the MRC-WPT system under the 
general setup with multiple transmitters and/or receivers. The system with two transmitters and a single receiver or 
a single transmitter and two receivers has been studied in lfT2l — l[T6l . while their analytical results cannot be applied 
for a system with more than two transmitters/receivers. Furthermore, to our best knowledge, there has been no work 
on rigorously establishing a mathematical framework to jointly design parameters in the multi-transmitter/receiver 
MRC-WPT system for its performance optimization. 

In this paper, as shown in Fig. Q] we consider a point-to-multipoint MRC-WPT system, where one transmitter 
connected to a stable energy source sends wireless power simultaneously to a set of distributed receivers, each of 
which is connected to a given load. We extend the results in lfl2ll - lfT6l to derive closed-form expressions of the 
transmit power drawn from the energy source and the power delivered to each load, in terms of various parameters in 
the system. Our results reveal a near-far fairness issue in the case of multiuser wireless power transmission, similar 
to its counterpart in wireless communication. Particularly, a receiver that is far away from the transmitter and thus 
has a small mutual inductance with the transmitter generally receives lower power as compared to a receiver that 
is close to the transmitter. We then show that the near-far issue can be optimally solved by jointly designing the 
receivers’ load resistances to control their received power levels, in contrast to the method of adjusting the transmit 
beamforming weights to control the received power in the far-held microwave transmission based WPT ifTTIl . tfT8l . 

Specifically, we first study the centralized optimization problem, where a central controller at the transmitter 
which has the full knowledge of all receivers, including their circuit parameters and load requirements, jointly 
designs the adjustable load resistances to minimize the total power consumed at the transmitter subject to the given 
minimum harvested power requirement of each load. Although the formulated problem is non-convex, we develop 
an efficient algorithm to solve it optimally. Then, for ease of practical implementation, we consider the scenario 
without any central controller and devise a distributed algorithm for adjusting the load resistances by individual 
receivers in an iterative manner. In the distributed algorithm, each receiver sets its load resistance independently 
based on its local information and a one-bit feedback shared by each of the other receivers, where the feedback of 
each receiver indicates whether the harvested power of its load exceeds the required level or not. Finally, through 
simulation results, it is shown that the distributed algorithm can achieve close-to-optimal performance as compared 
to the solution of the centralized optimization. 


II. System model 

We consider an MRC-WPT system with one transmitter and N receivers, indexed by n, n <E J\f = {!..... A r }, 
as shown in Fig. Q] The transmitter and receivers are equipped with electromagnetic (EM) coils for wireless power 
transfer. An embedded communication system is also assumed to enable information sharing among the transmitter 
and/or receivers. The transmitter is connected to a stable energy source supplying sinusoidal voltage over time 
given by v tx (t) = Re{u tx e J,i,t }, with v tx denoting a complex voltage which is assumed to be constant, and w > 0 
denoting the operating angular frequency of the system. Each receiver n is also connected to a given load (e.g. a 
battery charger), named load n, with resistance x n > 0. It is assumed that the transmitter and each receiver n are 
compensated by series capacitors with capacities c tx > 0 and c n > 0, respectively. Let i tx (t) = Rc{?' tx e ?, ' ;/ }, with 
complex i tx , denote the steady state current flowing through the transmitter. This current produces a time-varying 
magnetic flux in the transmitter’s EM coil, which passes through the receivers’ EM coils and induces time-varying 


3 



3D 

fD 

n 

fD 

< 

fD 


33 

fD 

O 

fD 

<' 

fD 


Fig. 1: A point-to-multipoint MRC-WPT system with communication and control. 


currents in them. We thus denote i n (t) = Rc{i n e^ wl }, with complex i n , as the steady state current at receiver n. 

We denote r tx > 0 ( r n > 0) and l tx > 0 (l n > 0) as the internal resistance and the self-inductance of the EM 
coil of the transmitter (receiver n), respectively. We also denote the mutual inductance between EM coils of the 
transmitter and each receiver n by h n > 0, with h n < where its actual value depends on the physical 

characteristics of the two EM coils, their locations, alignment or misalignment of their oriented axes with respect 
to each other, the environment magnetic permeability, etc. For example, the mutual inductance of two coaxial 
circular loops that lie in the parallel planes with separating distance of d meter is approximately proportional to 
d~ 3 0. Moreover, since the receivers usually employ smaller EM coils than that of the transmitter due to size 
limitations and they are also physically separated, we can safely ignore the mutual inductance between any pair of 
them. The equivalent electric circuit model of the considered MRC-WPT system is shown in Fig. [Q in which the 
natural angular frequencies of the transmitter and each receiver n are given by w tx = 1/ \fk x c lx and w n = 1/ \/l n c n , 
respectively. We set c tx = and c n = l ~ 1 w~ 2 , Vn G Af, to ensure that the transmitter and all receivers 

have the same natural frequency as the system’s operating frequency w, named resonant angular frequency, i.e., 
t^tx = w\ = ... = wn = w. 

We assume that the transmitter and all receivers are at fixed positions and the physical characteristics of their 
EM coils are known; thus, h n , Vn € AT, are modeled as given constants. We treat the receivers’ load resistances x n , 
Vn € AT, as design parameters, which can be adjusted in real-time lfl5l to control the performance of the MRC-WPT 
system based on the information shared among different nodes in the system via wireless communication. 

III. Performance Analysis 

In this section, we first present our analytical results. A numerical example is then provided to draw useful 
insights from the analysis. 


A. Analytical Results 

Define v = [t' tx - 0| X; .v] 7 and i = |i tx , % \,..., V\-], where v is the voltage vector and i is the current vector that 
can be obtained as a function of v. Let R = Diag(?q,..., r^v), X = Diag(xi,... , xn), and h = [h\, ..., Iin] T - 
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By applying Kirchhoff’s circuit laws to the electric circuit model in Fig. |T] we obtain 

T i _1 

—jwh 


i = 


v = A 1 v, 


r\x 

—jwh R + X 

where A G (£FV+i)xFV+i) is called the impedance matrix. The determinant of A is given by 


( 1 ) 


N 


N 


( 2 ) 


det(A) = (r tx + w 2 ^ h 2 k (r k + x k ) *) ( r k + x k ), 

k— 1 k= 1 

where it can be easily verified that det(A) > 0. Then, we define B = A” 1 , which is called the admittance matrix. 
Let B(b, l ) denote the element in row b and column l of B. We simplify dTJ as 


i = [B(l,l),...,B(N + l,l)) T v tx . 


(3) 


It can also be shown that B(b , 1), b G {1,... , N + 1}, is given by 

1 


B(b.l) = { 


r t x + w 2 T,k=i h U r k + Xk) 1 
wh b -i{r b i + x b -i)~ l 


/ r tx + w 2 J2k=i h U r k + x k )~ 1 

By substituting © into ©, it follows that 

1 


if 6 = 1, 
otherwise. 


(4) 


hx — 


r tx + w 2 J2k=i h K r k + x k ) 1 
wh n {r n + x n )~ l 


in =J- 


Vtx, 


Vtx, Vn G J\f. 


(5) 


( 6 ) 


'r t x + w 2 J2k=i h K r k + x k ) 1 

The power drawn from the energy source, denoted by p tx , and that delivered to each load n, denoted by p n , are 
then obtained as 

12 i 


Ptx = iRe{n, x i t ;}= ^--— 

2 2 r t x + w 2 Ylk=i h l( r k + x k )- 1 


1 ,• , 
Pn — -Xn\i n \ 


Ptx | 


w 2 h 2 n x n (r n + x r 


\—2 


(7) 


( 8 ) 


2 (’mWEiJLifclh + x*)- 1 ) 2 ’ 

where i* x is the conjugate of i tx • From it follows that the power delivered to each load n increases with the 
mutual inductance between EM coils of its receiver and the transmitter, i.e., h n . This can potentially cause a near-far 
fairness issue since a receiver that is far away from the transmitter in general has a small mutual inductance with 
the transmitter; thus, its received power is lower than a receiver that is close to the transmitter (with a larger mutual 
inductance). We accordingly define p sum = ^2^=1 Pn as the sum (aggregate) power delivered to all loads, where 
we always have p sum < p t x- 

In the following, we study impacts of changing the load resistance of one particular receiver n, i.e., x n , on the 
transmitter power p tx , its received power p n and that delivered to each of the other loads rn G N, m / n, i.e., p m , 
as well as the sum power delivered to all loads p sum , assuming that all other load resistances are fixed. 


Property 1. p tx strictly increases over x n > 0. 

This result can be explained as follows. From ©, it is observed that the transmitter current i tx strictly increases 
over x n > 0. Flence, due to the fact that the energy source voltage n tx is fixed, it follows that p tx given in ([7]) 
strictly increases over x n > 0. 
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Property 2. p m , Vm / n, strictly increases over x n > 0. However, p n first increases over 0 < x n < x n , and then 
decreases over x n > x n , where 

Xn = {r n (r tK + fi n ) + w 2 h 2 n )/(r ix + fi n ), (9) 

with fin = w 2 Ysk^n h l( r k + Xk)^ 1 - 

The above result can be justified as follows. From ©, it follows that for each receiver rn, rn f n, its current 
i m strictly increases over x n > 0. This is because i tx increases with x n , and as a result, i m increases due to the 
mutual coupling between EM coils of receiver m and the transmitter. Hence, the received power p rn defined in © 
also strictly increases over x n > 0. On the other hand, it follows from © that for receiver n, its current i n strictly 
decreases over x n > 0. Moreover, from (©, it follows that the decrement in \i n \ 2 is smaller than the increment of 
x n when 0 < x n < x n ; thus, p n increases with x n in this region, while the opposite is true when x n > x n . 

Property 3. If r tx + fi n — 2 <p n < 0, p sum strictly increases over x n > 0, where <p n = w 2 Jfk^n hl x k( r k + Xk)~ 2 ; 
otherwise, p sum first increases over 0 < x n < x n , and then decreases over x n > x n , where 

Xn = (fn(x tx T fin) T W h n T < 2v n l pn) / (ft\ T fin 2(fi n ) . (10) 

This property is a direct consequence of Property 2. 

B. Numerical Example 

We consider an MRC-WPT system with N = 3 receivers, where v tx = 25\/2V, r tx = 0.35H, f tx = 6.35//H, 
r n = 0.150, l n = 0.85/rH, Vn € A f, h = [2.3,1.1,0.9]/iH, and w = 2.2 x 10 6 rad/s. In this example, receiver 1 is 
closest to the transmitter and thus it has the largest mutual inductance, while receiver 3 is farthest. For the purpose 
of exposition, we fix X 2 = £3 = 7.50. We plot p tx , p n , Vn € N, and p sum , versus the resistance of load 1, x\, in 
Fig. |2l It is observed that p tx , P 2 , P:>, and p sum all increase over x,\ > 0. Note that in this example, the condition 
^ + fin~ 2 ip n < 0 holds in Property 3. However, p\ first increases over 0 < x\ < x\ = 15.80, and then declines 
over xi > 15.80. These results are consistent with our above analysis. Finally, we point out that changing x\ not 
only affects p\, but also the power delivered to other loads. For instance, receiver 1 can help receivers 2 and 3, 
which are farther away from the transmitter, to receive higher power by increasing x\. This is a useful mechanism 
that will be utilized to solve the near-far issue. 
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IV. Centralized Optimization 

In this section, we optimize the receivers’ load resistances x n , Vn € A f, to minimize the power drawn from 
the energy source at the transmitter subject to the given load constraints. We assume a central controller at the 
transmitter, which has full knowledge of the receivers, including their circuit parameters and load requirements, to 
implement the proposed centralized optimization. 


A. Problem Formulation 


We assume that the resistance of each load n can be adjusted over a given range x n < x n < x n , where x n > 0 
and x n > x n are lower and upper limits of x n due to practical considerations. We also assume that the power 
delivered to each load n should be higher than a certain power threshold p n > 0. Hence, we formulate the following 
optimization problem to minimize the power drawn from the energy source at the transmitter. 

kx| 2 i 


(PI): min _ 

{x n <x n <x n } 


s.t. 


Kp 


2 r tx + w 2 J2k=i h K r k + x k) 1 

w 2 h^x n (r n + x n y 2 


> p , Vn € J\f. 

— r\ ' 


2 {rtx + vJ 2 J2k=i h l( r k + Xk) X ) 2 

(PI) is a non-convex optimization problem. However, in the next we propose an efficient algorithm to solve (PI) 
optimally. 


B. Proposed Algorithm 

We define an auxiliary variable z = l/(r tx + w 2 ^k( r k + x k )~ l ) > 0. Since x n < x n < x n , Vn € A f, we 

have z < z < z, where z = l/(r tx + w 2 J2k=i h K r k + x k ) X ) and z = l/(r tx + w 2 J2k=i h K r k + x k ) : )- Then, 
we rewrite (PI) as 


(P2) : min |ftx| 2 t/2 

{linl^Xri^Xn}-, Z<Z<Z 


s.t. z 2 w 2 h 2 x n (r n + x n ) 2 >P, VneAC 

Z n 

(11) 

N 


Ax + w 2 ^ h 2 k (r k + Xk)- 1 = z- 1 . 

(12) 


k =1 


Although (P2) is still non-convex, we can solve it in an iterative manner by searching for the smallest z, z< z <~z, 
under which (P2) is feasible. Staring from z = z, we test the feasibility of (P2) given 2 by considering the following 
problem. 

(P3) : Find {x n < x n < x n , s.t. (Hill and (fl2]) |. 

If (P3) is feasible, then we set the optimal objective value of (P2) as z, which can be attained by any feasible 
solutions to (P3). Otherwise, we set z = z + Az, where A z > 0 is a small step size. We repeat the above procedure 
until (P3) becomes feasible or z >~z. The following proposition summarizes the feasibility conditions for (P3). 

Proposition 1. Given z, with z<z<~z, (P3) is feasible if and only if all conditions listed below hold at the same 
tune: 

Cl: z > 2y/r n /a n , Vn € A f, where a n = \v tx \ 2 w 2 h‘} l /(2p n ). 
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TABLE I: Algorithm for optimally solving (PI). 

Algorithm 1 

a) Given x n > 0 and x n > x n , Vn G Af, compute z and z. Initialize z z, A 2 > 0, and Flag «— 0. 

b) While z <~z and Flag = 0 do: 

1) Given z, check the conditions listed in Proposition |T| 

2) If at least one condition does not hold, then set z = z + A z. Otherwise, set Flag = 1. Choose any (y i,... ,y n ) G where the set <I> is 
given in condition C3 of Proposition |T| Set x n = 1 /y n — r n , Vn G A f. 

c) If Flag = 1, then return (xi ,... ,xn) as the optimal solution to (PI). Otherwise, problem (PI) is infeasible. 


C2: < x n < x^ and/or x\{ < x n < x^, Vn € Af, where xf/ = {a n z 2 /2 — r n ) — z\/ a n (a n 2: 2 /4 — r n ) and 

x v n = {a n z 2 /2 - r n ) + z^a n (a n z 2 /4 - r n ). 

C3: <4> = {(yi,... ,y N ) \ V n < Vn < y n , Vn € Af, r tx + w 2 h lvk = z~ 1 } + 0, where y_ n = 1 /{r n + 

min{x, x^}), and y n = 1 /{r n + max{x, x^}). 

Given any (yi,... ,y n ) € 't, where <I> is given in C3 of Proposition [Q the corresponding feasible solution to 
(P3) is obtained by a change of variable as x n = l/y n — r n , Vn € Af. Note that the obtained (x\,... ,xn) solves 
(PI) optimally. To summarize, the algorithm to solve (PI) is given in Table 1, denoted by Algorithm 1. 

V. Distributed Algorithm 

In this section, we present a distributed algorithm for (PI), where it is suitable for practical implementation 
when a central controller is not available in the system. In this algorithm, each receiver adjusts its load resistance 
independently according to its local information and a one-bit feedback from each of the other receivers indicating 
whether the corresponding load constraint is satisfied or not. We denote the feedback from each receiver n which 
is broadcast to all other receivers as FB n € {0,1}, where FB n = 1 ( FB n = 0) indicates that its load constraint 
is (not) satisfied. 

In Section [Dll we show that the power delivered to each load n, i.e., p n , has two properties that can be exploited 
to adjust x n . First, p n strictly increases over x rn > 0, Vm / n, which means that other receivers can help boost p n 
by increasing their load resistances. Second, p n has a single peak at x n = x n , assuming that other load resistances 
are all fixed. Thus, over 0 < x n < x n , receiver n can increase p n by increasing x n ; similarly, for x n > x n , it 
can increase p n by reducing x n . Although receiver n cannot compute x n from d9j) directly due to its incomplete 
information on other receivers, it can test whether 0 < x n < x n , x n = x n , or x n > x n as follows. Let p n (x+), 
Pn(x n ), and p n (x~) be the power received by load n when its resistance is set as x n + Ax, x n , and x n — Ax, 
respectively, where Ax > 0 is a small step size. Assuming all the other load resistances are fixed, receiver n can 
make the following decision: 

• If Pn(x+) > Pn (Xn) and Pn(x~) < Pn{x n ), then 0 < Xn <Xn, 

• If Pn(x+) < Pn{x n ) and p n (x~) < Pn(x n ), then x n = x n ll| 

• If Pn(x+) < Pn (xn) and Pn(x~) > Pn{x n ), then x n > x n . 

Next, we present the distributed algorithm in detail. The algorithm is implemented in an iterative manner, say, 

starting from receiver 1, where in each iteration, only one receiver n adjusts its load resistance, while all the other 
receivers just broadcast their individual one-bit feedback FB m , rn / n, at the beginning of each iteration. Initialize 

‘More precisely, if p n (x+) < Pn{x n ) and p n (x~) < Pn{x n ), then x„ - Ax < x„ < x„ + Ax. 







TABLE II: Distributed algorithm for (PI). 

Algorithm 2 

a) Initialize Itr = 1 and K max > 1. Each receiver n randomly chooses x n £ \x n , x n ]. 

b) Repeat from receiver n = 1 to n = N: 

— Receiver n receives FBm from all other receivers m ^ n and updates its load resistance x n according cases C1-C5. 

- If Itr = Ifmaxi then quit the loop and the algorithm terminates. 

- Set Itr = Itr + 1. 


by randomized x n € [x n . x n ], Vn € A f. At each iteration for receiver n, if p n < P n , then it will adjust x n to 
increase p n . To find the correct direction for the update, it needs to check for its current x n whether 0 < x n < x n , 
x n = x n , or x n > x n holds, using the method mentioned in the above. On the other hand, if p n > p n , receiver 

n can increase x n to help increase the power delivered to other loads when there exists any rn ^ n such that 

FB m = 0 is received; or it can decrease x n to help reduce the transmitter power when FB m = 1 , Vm / n. In 
summary, we design the following protocol (with five cases) for receiver n to update x n . 

Cl: If p n < p n and 0 < x n < x n , set x n <— min{x n , x n + Ax}. 

C2: If p n < p^ and x n > x n , set x n <— max{x n , x n — Ax}. 

C3: If p n > p , x n 7 ^ x n , and 3m / n, FB m = 0, set x n <— min{x n , x n + Ax}. 

C4: If p n > p , x n 7 ^ ± n , and FB m = 1, Vm / n, set x n <— max{x n , x n — Ax}. 

C5: Otherwise, no update occurs. 

In addition, we assume that there is a maximum number of iterations, denoted by K max > 1, after which 
the algorithm will terminate, regardless of whether it converges to a stable point (xi,..., xjv) or not. However, 
when the algorithm converges/terminates, the power constraints given in (ITTb may or may not hold for all loads, 
depending on the initial values of x n ’s. If constraint (fTTb holds for all loads, then the obtained (xi,..., xjv) is a 
suboptimal solution to (PI); otherwise, it is infeasible for (PI). The distributed algorithm is summarized in Table 
2, as Algorithm 2. 


VI. Simulation Results 

We consider the same system setup as that in Section IT1I-RI We set x n = 0.01H and x n = 10012, Vn € N. We 
also set p 1 = 250W, p 2 = 50W, and p, ; varying as 0 < p,, < 50W. Note that (PI) is feasible under the above 
setting. For Algorithm 1, we use As = 10~ 3 . For Algorithm 2, we use Ax = 10~ 3 and K m;iX = 10 5 , which is 
sufficiently large such that the algorithm converges to a stable point, while there is no guarantee that the power 
constraints given in (fill hold for all loads at this point. Therefore, to evaluate the performance of Algorithm 2, 
we averaged its result over 200 randomly generated initial points for each of which the algorithm converged to a 
feasible solution to (PI). In Fig. [3 we plot p tx versus /rj. It is observed that p lK obtained by Algorithm 1 is lower 
than that by Algorithm 2, while the gap is quite small, for all values of p^. This is expected since Algorithm 1 
solves (PI) optimally, while Algorithm 2 in general only returns a suboptimal solution. 

VII. Conclusion 

In this paper, we study a point-to-multipoint MRC-WPT system with distributed receivers. We derive closed- 
form expressions for the input and output power in terms of the system parameters. Similar to other multiuser 
wireless applications such as those in wireless communication and far-field microwave based WPT, a near-far 








fairness issue is revealed in our considered system. To tackle this problem, we propose a centralized algorithm 
for jointly optimizing the receivers’ load resistances to minimize the transmitter power subject to the given load 
constraints. For ease of practical implementation, we also devise a distributed algorithm for receivers to iteratively 
adjust their load resistances based on local information and one-bit feedback from each of the other receivers. We 
show by simulation that the distributed algorithm performs sufficiently close to the centralized algorithm with a 
finite number of iterations. As a concluding remark, MRC-WPT is a promising research area for which many tools 
from signal processing and optimization can be applied to devise new solutions, and we hope that this paper will 
open up an avenue for future work along this direction. 
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